What galaxy surveys really measure 
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In this paper we compute the quantity which is truly measured in a large galaxy survey. We 
take into account the effects coming from the fact that we actually observe galaxy redshifts and 
sky positions and not true spatial positions. Our calculations are done within linear perturbation 
theory for both the metric and the source velocities but they can be used for non-linear matter 
power spectra. We shall see that the complications due to the fact that we only observe on our 
background lightcone and that we do not truly know the distance of the observed galaxy, but only 
its redshift is not only an additional difficulty, but even more a new opportunity for future galaxy 
surveys. 
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I. INTRODUCTION 

All the photons which we receive have been emitted 
on our past light cone. In cosmology, looking far away 
always also means looking into the past. If the redshift 
of the objects under consideration is small, z <C 1, and 
evolution is relevant only on cosmological time scales, 
this effect is small and may be neglected. However for 
redshifts of order unity or larger, the fact that we are 
not observing a spatial hypcrsurfacc but a part of the 
background lightcone becomes relevant. 

If we observe the large scale distribution of galaxies, we 
usually compare the true, observed distribution with the 
one of an unperturbed universe with background density 
p and measure its fluctuations, <5(x) = (p(x)— p)/p, where 
p is usually the mean observed galaxy density. This is 
then cast in the power spectrum, 

(<5(k)<5(k')) = (2 7 r) 3 5(k-k')P (5 (fc) , 

where <5(k) is the Fourier transform of <5(x) and we as- 
sume statistical homogeneity and isotropy. For small 
galaxy catalogs one may assume that we measure the 
density fluctuation today, <5(x) = <5(x, to), but already 
for the Sloan Digital Sky Survey (SDSS) which deter- 
mines the galaxy distribution out to z ~ 0.2 or 0.5 (for 
Luminous Red Galaxies, LRG's) it is no longer a good 
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approximation, to compare the observed power spectrum 
with the above defined Pa (to). Time evolution of Pg can 
be taken into account by multiplying the power spectrum 
with a growth factor. In addition to this there is the is- 
sue of gauge. The density fluctuation £(x, t) which we 
calculate in a given Friedmann background is not gauge 
invariant. It depends on the background Friedmann uni- 
verse we compare the observed p(x, t) with. This is the 
cosmological gauge problem [l[ . 

There are several attempts in the literature to deal 
with these issues, but they are so far incomplete. People 
have considered individual observational effects like red- 
shift space distortions Q, the Alcock-Pacinski effect Q 
or lensing. A first full treatment is attempted in In 
the present paper we shall go beyond this work and de- 
termine the spectrum truly in terms of directly observ- 
able quantities. We derive gauge invariant expressions 
which are correct to first order in perturbation theory 
and which are straightforward to compare with observa- 
tions. This is an important first step for this problem as 
the gauge issue is mainly relevant on very large scales, 
where perturbations are small so that first order pertur- 
bation theory is justified. 

Our results will be most significant for future galaxy 
catalogs like BOSS @, DES 0], PanStarrs Q or, espe- 
cially Euclid [1], but also an analysis of SLOAN-7 @ 
along the lines outlined here is interesting. 

Notation: Wc work with a flat Friedmann back- 
ground and in conformal time, t, such that 

ds 2 = a 2 (t) (-dt 2 + (JydxV) . 
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A photon geodesic in this background which arrives at 
position xo at time to and which has been emitted at 
affine parameter A = at time ts, moving in direction n 
is then given by (x M (A)) = (A + ts, x + (A + t s - £o)n). 
Here A = t — ts — rs — r, where r is the comoving 
distance r = |x(A) — xo|, hence dr = —d\. We can of 
course choose xq = 0. 



II. THE MATTER FLUCTUATION SPECTRUM 
IN REDSHIFT SPACE 

In a galaxy rcdshift survey, the observers measure the 
number of galaxies in direction n at redshift z, let us 
call this iV(n, z)dtt n dz. They then average over angles 
to obtain their redshift distribution, (N)(z)dz. From this 
they can build directly the redshift density perturbation 
[22I ] i.e. the perturbation variable 

p(n,z) - (p){z) V (n,z) V{z) 



S z (n,z) = 



N(n,z)-{N)(z) 6V(n,z) 



(1) 



(N)(z) V(z) ■ 

Here V(n, z) is the physical survey volume density per 
redshift bin, per solid angle. The volume is also a per- 
turbed quantity since the solid angle of observation as 
well as the redshift bin are distorted between the source 
and the observer. Hence V(n, z ) = V{ z ) + 8V(n, z )- The 
truly observed quantity is the perturbation in the number 
density of galaxies 



N(n,z)-(N)(z) 



8 z (n,z) 



5V(n,z) 



A(n,z) (2) 



(N)(z) ^ '~" V(z) 

which therefore must be gauge invariant. Actually, as 
we shall see, both 5 z (n,z) and SV(n, z)/V(z) are gauge 
invariant. This is not surprising, as we could measure the 
volume perturbation also with other tracers than galaxies 
and it is therefore measurable by itself and hence gauge 
invariant 

We neglect biasing in our treatment as we want to 
keep the expressions as model independent as possible. 
We shall add only some comments on how simple linear 
biasing could be included. 



Here z = z(t) is the redshift of a background Fried- 
mann universe we compare our perturbation with and 5z 
is the redshift perturbation to this universe. Moreover 
p(n, z(t)) = p(t) + Sp(n,t), where the time is obtained 
by solving the background relation z = z(t). Note that 
p(z) = p(z + Sz) deviates to first order from p(z). Cleary, 
both 5z and 5p depend on the chosen background and 
are hence gauge dependent. However their combination 
in Eq. (J3J) must turn out to be gauge invariant, as it is in 
principle observable. 

Let us first compute the redshift in a perturbed Fried- 
man universe with metric 

ds 2 = a 2 {t)\- {l + 2A)dt 2 -2B i dtdx i + (4) 

+ [(1 + 2H L )hj + 2H Tij + -HI,, d.r'd.r' . 

Here Hij is the transverse traceless gravitational wave 
term and A, Bi, Hl and Hrij are scalar degrees of free- 
dom, two of which can be removed by gauge transfor- 
mations. In Fourier space Bi = —kiB and ifrij = 
(kikj — 8ij/3)HT- For simplicity, we shall neglect the 
contribution from gravitational waves in the main text. 
In the appendix we include also these terms. We consider 
a photon emitted from a galaxy, the source, S, which is 
moving in direction n (hence, to lowest order, it is seen 
under the direction — n from the observer O). We denote 
the peculiar velocities of the source and observer by vs 
and vo . The observer receives the photon rcdshifted by 
a factor 



1 + z 



(n ■ u)j 
(n ■ u) c 



(5) 



We solve the equation for the photon geodesic n = 
a~ 2 (l + 5n°, n + Sn), where n denotes the unperturbed 
photon direction at the observer. Using that u = a _1 (l — 
A, v), where v is the peculiar velocity, we find by the 
same calculation which leads to Eq. (2.228) in [1[ (see 
also [H) 



1 + z 



ajto) 
a(t s ) 



Hr 



\h t + n • V 



$ + * 



Is 



($ + *)cZA} 



(6) 



A. Computation of S z (n, z) 

Let us first relate 8 z (n,z) to the well known gauge 
dependent quantity 8(x, t). For this we note that to first 
order 



S z (n,z) 



p{n,z) - p(z) = p(z)+6p(n, z) - p{z) 

p{z) p(z) 
p(z — 8z) + Sp(n, z) — p(z) 

Sp(n,z) dpSz(n,z) 
p{z) dz p(z) 



(3) 



The first term is simply 1 + z. Here t denotes confor- 
mal time, $ and $ are the Bardccn potentials and V 
is the gauge invariant velocity perturbation which corre- 
sponds to the ordinary velocity perturbation in longitu- 
dinal gauge. For more details see [l|,El] an( i Appendix \M 
In this redshift perturbation the dipole term n-V(xo, to) 
is the only term in the square bracket in ^ which de- 
pends on directions when evaluated at xo. The terms 
evaluated at the emission point of course do all depend 
on n via the position of the emission point which, to 
lowest order, is simply X5 = xq — n(to — t[zs))- The in- 
tegral extends along the unperturbed photon trajectory 
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from the emission point, where we set A = 0, to our po- 
sition where A = to — t$ = rg. Eq. © implies that the 
redshift perturbation is 



Sz 



-(1 + z) (H L + -H T + n • V + $ + *) (n, z) 



($ + <&)dA 



(7) 



where we have neglected the unmeasurable monopole 
term and the dipole term from the observer position. We 
indicate the source position by the direction it is seen 
under, — n, and its observed redshift z. To lowest order 
x(n, z) = — rs(z)n. To obtain the density fluctuation in 
redshift space, we now use that ^| = 3j^. With this we 
obtain 

S z (n,z) = £> 9 (n,z) + 3(V-n)(n,z)+3(tt + $)(n,z) 



-3 / (V + <f>)(n,z(t))dt 

'ts 



(8) 



Here we relate a perturbation variable in direction n at 
redshift z to its unperturbed position and time, /(n, z) = 
f (x(n, z), t(z)) and overdots arc partial derivatives with 
respect to t, the second argument in /(x, t). 

D g is the density fluctuation on the uniform curvature 
hypersurface. It is related to the density fluctuation in 
co-moving gauge, D cm by [l|, [l(J 

Dcm = D = D g + 3$ + Zk^HV. 

If we would want to introduce a bias between the matter 
density and the galaxy density it would probably be most 
physical to assume that both galaxies and dark matter 
follow the same velocity field as they experience the same 
gravitational acceleration. We then expect that biasing 
should be applied to the density fluctuation in co-moving 
gauge, D cm , not to D g . On small scales such differences 
are irrelevant but on large scales they do become relevant 
as becomes clear when considering the (linear) power 
spectra for the different density fluctuations variables, 
see Fig. Q] 
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FIG. 1: The (linear) matter power spectrum on the uni- 
form curvature hypersurface (top curve, green), in longitudi- 
nal gauge (middle curve, red) and in co-moving gauge (bottom 
curve, blue). 



We want to express the volume element in terms of the 
polar angles at the observer position, 9o and ipo , and the 
observed redshift z. We have 



dV 



'9 i-nvafrW 



,8x v 8x a dx@ 



dz 80s dips 
v{z,0o,ipo)dzd9 o dipo , 



d(0o,<Po) 



dzd^odi^o 
(10) 



where we have introduced the density v which determines 
the volume perturbation, 



SV 
V 



5v 
v 



is the determinant of the Jacobian of the trans- 



d(6 s ,<ps) 
d{8o,<Po) 

formation from the angles at the source to the angles at 
the observer. Eq. (fTU)) is still exact. In a homogeneous 
and isotropic universe geodesies are straight lines and 
&S = 6o an d ips = <Po- In a perturbed universe the 
angles at the source are perturbed with respect to the 
angles at the observer and we have 0$ = 0o + 89 and 
= fo + Sip. Hence to first order the Jacobian deter- 
minant becomes 



d(0 Sl Ps) 



d{0 o ,<Po) 



1 850 ddp 



B. Volume perturbations 

As next step we compute the volume perturbation 
5V/V in Eq. which must be gauge invariant since 
also 5 Z is gauge invariant by itself. This is not surprising 
as it would in principle be a measurable quantity if we 
would have an 'unbiased tracer' of the volume. 

We consider a small volume element at the source po- 
sition. By this we mean the spatial volume seen by a 
source with 4- velocity it M . This is given by 



dV 



'9 e nva/3 U 



>*dx u dx a dxP. 



(9) 



Using the expression for the metric determinant, \J—g = 
a 4 (l + A + 3Hl) and the 4- velocity of the source, u = 
i(l — A,v l ), we find to first order 



a 3 (l + A + 3H L ) 



dr , . „ / 859 85p 



dr dt 
A— + v r — I r" sm« 
dz dz 



(12) 



Here dr / dz is to be understood as the change in comoving 
distance r with redshift along the photon geodesic. At 
linear order we can write (the distinction between z and 
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z is only relevant for background quantities) 

dr df dSr ddz df ( df dSr dSz df\ dt 
dz dz dz dz dz \ dt dX dX dz J dz' 

(13) 

where we have used that for first order quantities we can 
set dt = dX when we have to take the derivative along 
the photon geodesic. The last term of Eq. (fT3f contains 
the redshift space distortion which will turn out to be the 
biggest correction to the power spectrum. To lowest or- 
der along a photon geodesic —df/dz = dt/dz = —H^ 1 = 
—a/H, where H is the physical Hubble parameter and 
H = aH is the comoving Hubble parameter. With this 
the volume clement becomes 



a 4 r 2 sin 9q 



n 



d 

I cot 9 Q + — 



28r 
r 



dSr 
~dX 



a dSz 
U~dX 



dSip 
dip 

(14) 



To obtain the fluctuation of v wc subtract the unper- 
turbed part v(z). Note, however that we evaluate this at 
the observed redshift, z = z + Sz. Hence 

v(zj = v(z) + ^poz. 

From the unperturbed expression with a = l/(z + 1), 

?2 



v(z) = 



sin uor' 

[i + zyn 



one infers 



dv ._. / A 2 H \ 1 
— =v(z) [-4+ — + — 2 



dz 



f s n n 2 i + z 



(15) 



(16) 



Combining Eq. (|T4|) and ([TBI we obtain for the pertur- 
bation of the volume element 



Sv viz) — v(z) 

— (n,z) = -— = 

V v(z) 

( n d \ _ „ dScp 26r 
3H L + [cot0 + - 56+-^--vn+ — 



(17) 



dSr 



1 dSz 



H \ Sz 



dX H(l + z) dX 



ru n 2 i + z 



In order to express these quantities in terms of the 
perturbed metric and the peculiar velocity of observer 
and emitter, we need to compute the deviation vector 
that relates the perturbed geodesic to the unperturbed 
one &r M (A) = x^(X) — x^(X). Here we give only the main 
steps. More details on the derivation can be found in the 
appendix. We use 



~df 



dx* 1 dX 
~~d~X~dt 



(18) 



which leads to 

x°(ts) — —(to—ts)=rs at every order (19) 

x'fa) = -(t -*s)n*- / dX(5n l - tfSn ) (20) 

Jo 

to first order. 

In the following we neglect perturbations at the ob- 
server position since, as already mentioned, they give rise 
only to unmeasurable monopole term or a dipole term. 
Using the null geodesic equation for n M wc find 

Sx^ts) = + I dx(h m n a + h 0a n [ n a ^ (21) 



'o 





where r(A) = A. From this we obtain 

2 



Sr 



-8x l rij 



dXh a pn a n 13 



rs Q 



1 fdH T 



k k 2 \ dX 



2Hr, 



(22) 



dX dt 



I and r s = t Q - t s 



We also use that n l di + d t 

to lowest order. For the derivative of 5r we obtain 

d5r , . 1 dB 1 (d 2 H T dHr 

= -($ + *) H 1 - 2 - 

dX ' ' k dX k 2 \ dX 2 dX 

Similarly we find for the perturbed angles 
8x l e ftl 1 



(23) 



SO 



5(p 



rs 



rs Jo 



dX 



x [h ai n a e l e + -(rs - r)h a p. i e , e n a n f 



1i e 



1 



dX 



r s smV r s smV Jo 

XI 



(h m n a e l ip + -(r s - r)h a p^e) p n a 'i 



(24) 



(25) 



We have used that n l eei = n*e V i = 0. The second term 
of the integral in Eq. (f2~4"|) can be rewritten as 

1 



h a pje l g n a nP = -de(h a f})n a n^ 

= - de(h ali n a h si ) - h a ^de(n a n fi ) 



(26) 



where don a = —e l g 6i a , and analogously for ip. The angu- 
lar contribution to the volume then reads 

(cot 9 + d g )56 + djp = [ S dX - ^ 
Jo 2r sr 

1 



x 



cot Ode + dl 



rs 1 

dX- 



r 

.\ r r s 



sin' 

(cot 9 + de)h ia e l e n 



h a pn a n^ 




sm ( 



dX {rs r) A fj ($ + ^) 
r 

A n H T {t s ) 



k 2 r% 



(27) 
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where Aq denotes the angular part of the Laplacian 



A n = cot 00 fl + 02 



sin 2 0^ 



(28) 



It is interesting to note that the angular part of the vol- 
ume perturbation is not a gauge-invariant quantity by 
itself. If Ht +~ the angular and radial directions are 
mixed in a non-trivial way. This is not really surprising 
since the angular volume distortion is not a measurable 
quantity by itself. On the other hand, the convergence 
k (or the magnification fi) that are observable, contain 
in addition to the angular volume distortion other per- 
turbations (see [Hj]) and are consequently gauge in- 
variant. 

The redshift contribution to the volume perturbation 
is obtained by differentiating Eq. ([7}. 



f 



dSz 



H(l + z) d\ 
$ + * + H L 
1 



(29) 



Vn 



dA($ + *) 



3 jo 

dH L 1 dH T d(V ■ n) 



3 dX 



dX 



Putting everything together we find after several integra- 
tions by part and a total Laplacian of Ht which cancels 
a factor 1/k 2 , the following expression for the volume 
density perturbation 



d r V- 



dA($ 




d(V • n) 



dX 



(30) 



Here and in the following, the functions without argu- 
ment are to be evaluated at the source position xs = 
xq — n(to — ts) and at the source time ts- More details 
of the derivation of this result are given in the appendix. 

Adding the results © and (|5D)| we obtain the galaxy 
number density fluctuation in redshift space as defined 
in Eq. © 



A(n,z) 




* + <9 r (V • n) 



dA($ + r 



($ + *). (31) 



Here we have used that also pressureless matter moves 
along geodesies so that 

n- V + Hn- V-d r # = 0. 



Equation (j3"Tj) together with © and (|3U|) are our first 
main result . 

The first term in (|3ip is the gauge invariant density 
fluctuation. D g is the density fluctuation in the flat slic- 
ing. It is related to the density perturbation in Newto- 
nian gauge by D g = D s — 3$. In terms of D s the first 
three contributions combine to Dg + Q+ty = D a — 2$+\I r . 

The term H~ 1 d r (n- V) is the redshift space distortion. 
As we shall see in the next section, this is the largest 
single correction on intermediate scales. The second line 
comes from the redshift perturbation of the volume. It 
contains a Doppler term and the ordinary and integrated 
Sachs- Wolfe terms. The third line represents the radial 
and angular volume distortions. The second term in the 
integral on the third line is especially relevant on large 
scales, it is the lensing distortion. 



III. THE ANGULAR POWER SPECTRUM OF 
THE GALAXY DENSITY FLUCTUATIONS 

For fixed redshift A(z,n) is a function on the sphere 
and it is most natural to expand it in spherical harmonics. 
Let us do this with the result (|3"TT) 



A(n,z) =Y, a im(z)Y £m (n), C e (z) = (\a em \ 2 ). (32) 



Cm 



The coefficients ag m (z) are given by 

ai m {z) = / dO n F/ m (n)A(n, z). 



(33) 



The star indicates complex conjugation. 

The different terms in A(n, z) are either a perturbation 
variable evaluated at the source position or an integral 
of a perturbation variable over the unperturbed photon 
trajectory. Let us first consider a contribution from a 
perturbation variable at the source position, e.g. VP. We 
want to relate the Cg(z) spectra to the usual power spec- 
trum P^(k,t) which is defined by 

(*(k,t)**(k',t)> = (27r) 3 5(k-k')P*(fc,i). 

The delta function and the fact that depends only on 
the modulus of k, k = |k|, are a consequence of statistical 
homogeneity and isotropy. Expressing ^ in terms of its 
Fourier transform, 

*(x,t) = ^j J d 3 k*(k,t)e-^ *\ 

a short calculation (see e.g. [l[) gives 

■i r 

a f™( z s) = 2^2 J d 3 fcj£(fcr s )*(M s )y/ m (k). (34) 

Here jj is the spherical Bessel function of order £, 
see [13j . Correspondingly the contribution from an in- 
tegral J Q rs /(x(A), t(X))dX becomes 

4lizs) = ^[ S ^ J d 3 k Je (kX)f(k,t)Y; m (k). (35) 
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For a velocity term V • n we use that V(k, t) = ikV, so 
that V • nexp[i(k • n)r] = Vdkr exp[i(k • n)r]. With this 
one obtains 



= ^3 \ d 3 kj>(kr s )V(k,t)Y; m (k). (36) 



Vn/ 

Hm I 



The prime in jg denotes derivation w.r.t the argument. 
Finally, for the rcdshift space distortion, d r (V ■ n) 

ii ■ V(V • n) we have to use the above identity twice 
and arrive at 



a r (Vn) 



(zs) = d A kj'l{kr s )k- L V{k,t)Yl m {k). 

(37) 

One can now write down the C^(z)'s for one's theory 
of choice for the background and the perturbations, e.g. 
for modified gravity or a quintessence model. 

So far the derivation has been completely general. We 
have not used Einstein's equation. The only assumptions 
are that galaxies follow the distribution of matter which 
is made out of non-relativistic particles which move along 
geodesies, and that photons move along null geodesies. 

To proceed further, we have to be more specific. Here 
we just study the simplest model of purely scalar adia- 
batic perturbations, which have been generated at some 
early time in the past (e.g. inflation). If there are more 
e.g. isocurvature modes present, the subsequent calcula- 
tion has to be repeated for them. 

In the case of one adiabatic mode, all the perturbation 
variables arc given by transfer functions from some initial 



random variable that we take to be the Bardccn potential 
, 5. Hence 



*(M) 
d>(k,t) 

D g (k 7 t) 
V(k,t) 



T*(M)*in(k) 

T*(M)*ia(k) 
7b(M)*in(k) 

T v (M)* in (k). 



(38) 
(39) 
(40) 
(41) 



The transfer functions T. depend on the matter content 
and the evolution history of the Universe and on the the- 
ory of gravity which relates matter and metric degrees 
of freedom. What is important for us is that they are 
deterministic functions and do not depend on directions 
of k. We characterize the initial power spectrum by a 
spectral index, n, and an amplitude, A, 



fc 3 (* in (k)**„(k')) = {2nf5{k-k!)A{kt ) n - 



(42) 



We have introduced present time, to, in order to keep A 
dimensionless. From the CMB observations we know that 
it is of the order of A ~ 10 -8 . With these identifications 
we can now relate Ci{z) to the initial power spectrum 
A(kto) n ^ 1 ■ Inserting the above in expression (|3"Tj) a short 
calculation gives 



zs) 



(43) 



with 



F e (k,zs) = ji(kr s ) 



S j e (kX) (2 + ^— W + 1)1 (7* + T<s,)d\ + 
rs Jo V A / 



je(kX)(Ty + T<s,)d\. 




(44) 



Here rs = to — ts is the source position . 

We now evaluate and compare the amplitude of dif- 
ferent terms in a ACDM universe. Rather than entering 
in a precise numerical evaluation, we estimate the terms 
by using approximations for the transfer functions. This 
will help us to gain insight in the importance of the dif- 
ferent terms. We plan to do a full numerical evaluation 
which can be used to estimate cosmological parameters 
in future work. 

From the first order Einstein equations, neglecting 
anisotropic stresses from neutrinos, we can relate the 
transfer functions Tjj , Ty and 7$ to ■ We find 



Using the notation of [lj] (see also |15j). we decompose 
the transfer function T<s,(k,t) into a growth rate D\(a) 
and a time independent transfer function T(k) such that 



T*(M) = i^T(fc), 



(48) 



and we use CAMBcode [16J to compute T(k). The am- 
plitude of the power spectrum can be expressed as (l4l | 

A = r ^-5l ( p, . We choose Q m = 0.24, A = 



0.76 and erg 



D 1 (a=l) / 

0.75 leading to Sh 



5.6 • 10- 



2a 
2a k 



k 

Ho 



(45) 
(46) 

(47) 



A. The transversal power spectrum 

Let us first determine the C^'s at fixed rcdshift. They 
provide the transversal power spectrum, i.e. correlations 
on the sphere normal to the observer directions. Of 
course for the intrinsic density fluctuations these are not 
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FIG. 2: Top panel: The transversal power spectrum at (from 
top to bottom) zs = 0.1, zs = 0.5, zs = 1 and zs = 3. 
Bottom panel: The ratio between the new contributions (lens- 
ing+potential) and the total angular power spectrum at (from 
top to bottom) zs = 3, zs = 1, zs = 0.5 and zs = 0.1. Solid 
lines denote positive contributions whereas dashed lines de- 
note negative contributions. 

different from correlations in any other direction, but ob- 
servational effects on them are different. E.g., since we 
can only observe on the background lightcone, we can 
only see fluctuations on this sphere at the same time but 
not fluctuations which have a different radial distance 
from us. On the other hand, in general the same redshift 
does of course not imply the same look-back time, since 
both these quantities are perturbed in different ways. 

In Fig. [2] (top panel) we show the total transversal 
power spectrum at redshifts zs = 0.1, 0.5, 1 and 3. Note 
that the amplitude of the linear power spectrum from 
zs = 0.1 to zs = 0.5 is reduced by a factor 6 at t ~ 100 
and by a factor 20 at I ~ 10. This comes from the follow- 
ing fact: the transversal power spectrum is dominated by 
the density fluctuation and the redshift space distortion 
which are proportional to integrals of the form 

/f(£) V »>- 

At x = krs — £, this term goes like £ 4 and it is therefore 
expected to dominate at large I. However, since for a 
constant transfer function, this integral would diverge, it 
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FIG. 3: The dominant terms at redshifts (from top to bot- 
tom) = 0.1, 0.5, 1 and 3: density (red) , redshift space dis- 
tortion (green), the correlation of density with redshift space 
distortion (blue), lensing (magenta), Doppler (cyan), see Ta- 
ble U The potential terms are too small to appear on our 
log-plot. 
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FIG. 4: Top panel: The various terms as a function of zs for 
fixed value of t = 20: density (red), redshift space distortion 
(green), the correlation of density with redshift space distor- 
tion (blue), lensing (magenta), Doppler (cyan) and potential 
(black), see Table [U Solid lines denote positive contributions 
whereas dashed lines denote negative contributions. 
Bottom panel: The ratio between the new contributions (lens- 
ing+potential) and the total angular power spectrum as a 
function of zs for fixed value of £ = 20. 



is dominated by the maximum of the transfer function 
which is roughly at fc Gq . Since for z ~ 0.5, k cq rs ~ 
£, ji{k cq rs) oc l/(fc oq r,s) 2 which therefore decreases like 
l/r|. Already this simple observation tells us that the 
amplitude of the transversal power spectrum at different 
redshifts might offer a possibility to constrain r$(z) and 
the growth factor, which both depend on cosmological 
parameters in different ways. On the other hand, this is 
complicated by non-linear effects and biasing which are 
not accounted for in this work. 

The different contributions to the power spectrum at 
different redshifts are shown in more detail in Fig. [3J For 
Zs = 1, we show the spectrum up to £ = 600 while for 
the other redshifts we stop at £ = 100 beyond which the 
structure does not change anymore. We denote by D the 
density term in co-moving gauge, 



D = D B 



3$ 



3^V. 
k 



by z the redshift space distortion, by L the lensing term, 
by V the Doppler terms and by 'F the gravitational po- 
tential terms (see Table fl] for a definition of each term). 



Density 
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red 


redshift space 






distortion 


K dr(V-n) 


green 


lensing 
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magenta 


correlation 
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blue 


Doppler 
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cyan 


potential 
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black 



TABLE I: The color coding the different terms of Eq. (J3TJ in 
the angular power spectrum of A(n, z) as shown in Figs. [3] 
to [5] and [7] [9] 1101 In addition to the term given in the sec- 
ond column, all its correlations with the terms in the lines 
above are also included. Only the most dominant correlation 
between density and redshift space distortion is shown sepa- 
rately in blue. In Figs. l5l and 1 101 the 'standard terms', i.e. the 
top three lines, are represented together as the blue line. 



Cj? D represents for example the contribution from the 
density term alone and Cf z the correlation between the 
density and redshift space distortion. Except from the 
correlation between the density and redshift space distor- 
tion that we represent individually, we include the corre- 
lations with the smaller contribution. Note that usually 
the correlations between lensing, Doppler and gravita- 
tional potential are negligible and, except when explic- 
itly specified we do neglect them. Therefore, when we 
plot for example the lensing term (magenta), it contains 
C\ L + 2C\ D + 2Cf z . The formulae for the dominant 
CiS are given in Appendix [Bj 

The lensing term scales like £^ and is in principle of the 
same order as the density and redshift space distortion 
terms. However it is given by an integral of the form (see 
Appendix IB1 



£ 2 (£+l) 2 



dX 



rs 



-je(kr) 



which does converge when integrated over k. It is there- 
fore dominated at k = £/r. (We have used Limber's 
approximation [l7j to evaluate this integral which we 
have tested numerically and found to be of excellent ac- 
curacy.) The contribution of the lensing term becomes 
more important at larger source redshift for small £. But 
it always remains subdominant in the transversal power 
spectrum. In the bottom panel of Fig. [2] we plot the ra- 
tio between the new contributions, i.e lensing term plus 
potential term, and the total angular power spectrum. 
We see that neglecting the new contributions for z$ < 1 
represents an error of no more than 0.1 percent, whereas 
for zs = 3 the error amounts to a few percent. Note that 
we do not include the Doppler terms in the new contri- 
butions since they appear already in the original Kaiser 
formula (i~8l | (even though there the term from expan- 
sion oc T-t/Ti 2 . which is of the same order for redshifts 
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FIG. 5: Different terms of Ce(zs, at I — 20 for redshifts 
(from top to bottom) zs — 0.1, 0.5, 1 and 3, plotted as 
a function of z$> ■ standard term, i.e. Cf D + C| z + 2C<P Z 
(blue), lensing (magenta), Doppler (cyan), potential (black), 
see Tabled] Solid lines denote positive contributions whereas 
dashed lines denote negative contributions. 
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FIG. 6: The effect of a window function on the density 
contribution (top panel), redshift space distortion (middle 
panel) and lensing contribution (bottom panel). We have 
chosen zs = 0.1 and Azs = (no window, top curve, red), 
Azs = 0.002 (middle curve, green) and Azs = 0.01 (bottom 
curve, blue). 

z > 1 is not considered). In Fig. [4] (top panel) we de- 
pict the redshift dependence of all the terms for a fixed 
value of I = 20. The lensing and potential terms are both 
negative at small redshift and become positive at large 
redshift. This is due to the fact that at small redshift the 
dominant contribution comes from their correlation with 
the density that is negative, whereas at large redshift the 
dominant contribution is their auto-correlation, , re- 
spectively C**. The bottom panel of Fig. 0] shows the 
ratio between the new contributions and the total angu- 
lar power spectrum. The error induced by neglecting the 
new terms increases with redshift and it reaches a few 
percent at high redshift. 
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B. The radial power spectrum 
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FIG. 7: The effect of a window function with width Azs = 
O.lzs on the power spectrum d(zs) for redshifts (from top 
to bottom) zs = 0.1, 0.5, 1 and 3. The different curves are: 
density (red), redshift space distortion (green), the correlation 
of density with redshift space distortion (blue), lensing (ma- 
genta), Doppler (cyan) and gravitational potential (black), 
see Table |T] Solid lines denote positive contributions whereas 
dashed lines denote negative contributions. 



The results above give us the transversal power spec- 
trum at fixed redshift. But of course there is also a radial 
power spectrum which correlates fluctuation at different 
distances from us. This encodes different information 
and it is important to study them both. From the fact 
that the transfer function is not direction dependent, we 
infer that 

{atm(zs)a>t'm'(zs')) = Se,e'5m,m'Ce(zs,z S ') ■ (49) 
Hence the radial power spectrum is given by 
9 A C rlk 

C e (z s ,z s ,) = —J T (kto) n - 1 F e (k, Zs )F;(k,z s ,). 

(50) 

Here a interesting new phenomenon occurs: due to the 
fact that we evaluate F^{k,zs) at different redshifts, we 
also evaluate the Bessel functions je(krs) at different dis- 
tances r\g. This leads to a suppression of the result due 
to oscillations, if the region in /c-spacc where the inte- 
grand dominates has kr$ > £■ As we discussed above, 
this is the case for the fc 2 -term of the density fluctu- 
ations and for the redshift space distortion, the terms 
which dominate the transversal power spectrum. These 
terms are therefore substantially suppressed in the radial 
power spectrum. All other terms have convergent inte- 
grals of jj{krs), already when neglecting the turnover of 
the transfer function, hence they are suppressed by pow- 
ers of £ with respect to the lensing term. Therefore the 
lensing term dominates the radial power spectrum at low 
I. This is precisely what one sees in Fig. [5J where the 
lensing term (magenta) dominates for zs> significantly 
larger than zs- As in Fig. |4j at small redshift z$ = 0.1 
and zs = 0.5 the correlation density-lensing dominates 
(and is negative), whereas at large redshift zs = 3 the 
lcnsing-lensing term dominates. It is interesting to note 
how constant the lensing term remains while the density 
term and the redshift space distortion decay very rapidly 
with growing redshift difference. At zs = 1 the lcnsing- 
lensing term and its correlation with the density are of 
the same order of magnitude which explains the change 
of sign as zs> increases. Finally, at Zs = 0.1 the Doppler 
term dominates over the standard term for some very 
specific values of zs> ■ The first of them is actually the 
zero in the real space correlation function which e.g. at 
redshift = 0.1 corresponds to 6z = 0.011. 

An alternative way to measure radial correlations is to 
introduce a window function W(z, z') which corresponds 
to a smearing of fluctuations on scales smaller than some 
width Azs- We use a Gaussian window around some 
mean redshift zs with width Azs- This suppresses power 
which comes from values of k with kArs > t where 
Ars = r(zs + Azs) — r(zs)- This is also a more re- 
alistic case since we can measure the galaxy distribution 
only in redshift bins of some finite width. Already a small 
width does substantially affect the resulting spectrum of 
the density, see Fig. [H] top panel, and the redshift space 
distortion (middle panel). As expected the lensing term 
is insensitive to this smearing (bottom panel). 
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FIG. 8: Top panel: The total power spectrum at redshifts 
(from top to bottom) zs = 0.1, zs = 0.5, Zs = 1 and zs = 3 
smeared by a window function with width Azs = O.lzs. 
Bottom panel: The ratio between the new contributions (lens- 
ing+potential) and the total angular power spectrum at (from 
top to bottom) zs — 3, zs = 1, zs = 0.5 and zs = 0.1. Solid 
lines denote positive contributions whereas dashed lines de- 
note negative contributions. 

In Fig. [7] we show the effect of a 10% window on the 
different terms at different redshifts. As before, the terms 
which we indicate by 'lensing term' 'Doppler term' and 
'gravitational potential contributions' in the figure are 
not only the corresponding term themselves but also their 
correlations with all other terms. If the latter dominate 
such a contribution can become negative. For exam- 
ple the lensing contribution for zs = 1 changes sign at 
£ = 28. For £ > 28 it is dominated by negative cor- 
relations with the density while for £ < 28 the positive 
autocorrelation dominates. Since the power from scales 
smaller than fcArg is removed, the power at I truly cor- 
responds to that at k = £/r(z) in the power spectrum. 
The 'wiggles' in the density and in the velocity terms for 
zs = 0.1 are the baryon acoustic oscillations (BAOs), 
the first of which appears at £ ~ 15 for z$ = 0.1. They 
are also visible in the anti-correlation of the lensing term 
with the density for zs =0.5 and z$ — 1, but these terms 
are probably too small to be detected in real data. 

In Fig. [8] (top panel) we show the total Ce(zs)'s 
smeared with a 10% window function. Comparing it with 
Fig. [5J we mainly notice that the power is reduced sig- 
nificantly, by nearly 1.5 orders of magnitude. Furthcr- 



FIG. 9: The various terms as a function of zs for fixed value 
of I — 20 and smeared by a window function with width 
Azs = O.lzs: density (red), redshift space distortion (green), 
the correlation of density with redshift space distortion (blue), 
lensing (magenta), Doppler (cyan) and potential (black), see 
TablelU Here the correlations between the lensing and Doppler 
and the lensing and potential are not negligible at large zs, 
and they are included in the Doppler (cyan), respectively po- 
tential (black) curves. Solid lines denote positive contribu- 
tions, dashed lines denote negative contributions. 
Bottom panel: The ratio between the new contributions 
(lensing+potential) and the total angular power spectrum, 
smeared by a window function with width Azs = O.lzs, plot- 
ted as a function of zs for fixed value of £ = 20. 



more, at z$ = 0.1, the BAO's are clearly visible. In the 
presence of a window, different terms can dominate at 
different redshift and for different values of £. In the bot- 
tom panel of Fig. [5J we depict the ratio between the new 
contributions and the total angular power spectrum. Ne- 
glecting the new contributions induces an error of a few 
percent already at redshift 1, and this error increases to 
roughly 50 percent at redshift 3. Note that this ratio de- 
pends strongly on the width of the window function, and 
that a wider window would lead to a larger error. 

In Fig. IH1 (top panel) we plot the different terms as a 
function of redshift, for fixed value of £ = 20. Contrary 
to Fig. 31 where the lensing term always remains sub- 
dominant with respect to the density and redshift space 
distortion term, we see in Fig. [5] that for zs > 2.4 the 
lensing term dominates over the standard contribution. 
The redshift at which this dominance takes place depends 
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of course on the chosen window function: for larger Azs, 
the lcnsing term starts to dominate at smaller rcdshift. 
In the bottom panel of Fig. [9] we show the ratio between 
the new contributions and the total angular power spec- 
trum as a function of redshift for I = 20. From this figure 
we understand why in Fig. [5] (bottom panel), the ratio 
at zs = 1 is not significantly larger than at zs = 0.5. 
The lcnsing contribution changes sign around z$ = 0.9 
and consequently it is still small at zs = 1. At a rcdshift 
of zs = 1.5, however the error induced by neglecting the 
new terms is already of the order of 10 percent. 

In Fig. 1101 we show correlations between different red- 
shifts bins (with a 10% window function), for a fixed 
value of t = 20. As in Fig. [5] we see that the lensing term 
becomes dominant when the redshift separation between 
the bins increases. At large redshift, zs = 3, the lensing 
term is always dominant. The individual behaviour of 
each contribution is however quite different from Fig. [5] 
which is due to the smearing introduced by the win- 
dow function. Note that comparing the second panel in 
Fig. [10] with the results in [19| , we see that the redshift 
separation between their 4 different bins (their Fig. 13) 
is too small for the lensing contribution to be relevant. 
However, a similar measurement with one of the bins sit- 
uated around zs = 0.7 would already allow to detect the 
lensing contribution. 

Finally, we plot in Fig. [TT] the angular power spectrum 
integrated from the observer until a maximum redshift 
Zmax- This corresponds to the situation where the red- 
shifts of individual galaxies is unknown but obeys a given 
redshift distribution. Consequently only the integrated 
spectrum can be measured. We assume a flat distribution 
of galaxies between z m - ln = 0.1 and z max = 2 with Gaus- 
sian tails at both ends. In Fig. [TTJ we see that the only 
relevant contributions are the density and the lensing, 
more precisely the cross-correlation between the lcnsing 
and the density which is negative. The redshift space 
distortion contribution is, as expected, completely negli- 
gible when the galaxy redshifts are unknown. The lensing 
contribution, however is very relevant; it reduces the re- 
sult by roughly 40% of the contribution from the density 
alone. 



IV. CONCLUSIONS 

In this paper we have derived expressions for the 
transversal and radial galaxy power spectra, Ce(zs) and 
Ci(zs, zs 1 ), taking into account not only redshift space 
distortions, which have also been studied e.g. in [l9| but 
also all other relativistic effects to first order in pertur- 
bation theory. Within our accuracy we are in reasonable 
agreement with the simulated results of Ref. [lj| (their 
Fig. 4) which analyzes the SDSS data taking into account 
rcdshift space distortion but not the other terms, e.g. the 
lcnsing, appearing in our formula (|31[) . They also take 
into account non-linearities in the matter power spec- 
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FIG. 10: Cross-correlations between different redshift bins 
Ce(zs , Zs>) at £ = 20 with a 10% window function and plot- 
ted as a function of z' s . From top to bottom zs = 0.1, 0.5, 1 
and 3. Standard term, i.e. Cf D + Cf z + 2Gf z (blue), lens- 
ing (magenta), Doppler (cyan), potential (black), see Table [TJ 
Solid lines denote positive contribution whereas dashed lines 
denote negative contributions. 
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FIG. 11: Integrated power spectrum with a flat distribution 
between z m i n = 0.1 and z max = 2 with Gaussian tails at 
both ends. The density term is plotted in red and the lensing 
term in magenta. Note that the lensing term is completely 
dominated by its anti-correlation with the density and hence 
is negative. 



trum by using halofit [16j . This enhances their results 
with respect to ours. 

We have seen that by measuring Ci(z$, zs>) for differ- 
ent redshift differences and different £'s we can measure 
different combinations of terms which depend on cosmo- 
logical parameters in a variety of ways. Otherwise, one 
may measure the C/s smeared over a given redshift bin, 
Az s , 

C e (z s ,Az s ) = I dzdz'W{z,z')C e {z,z') 



where W is a window function centred at z$ with width 
Azs- Without smearing, the density contribution and 
the redshift space distortion always dominate. When 
smearing is included these terms are reduced and the 
lensing term can dominate. 

The method outlined in this paper represents a very 
flexible new path to estimate cosmological parameters 
and to test the consistency of the concordance model of 
cosmology. Of course, to do this we must master pos- 
sible degeneracies not only with biasing but also evolu- 
tionary effects which have not been discussed in this work 
and which may become relevant at redshift larger than 1, 
see [2(J for a discussion. A detailed parameter estimation 
forecast e.g. for Euclid is left as a future project. 
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Appendix A: Some details of the derivations for 

A(n,z) 

We consider a perturbed Friedmann metric, 
ds 2 = a 2 (t) (- (1 + 2 A)dt 2 -2B i dtdx i + (Al) 



[{l + 2H L )5 lJ +2H 



Tij 



2H l] ]dx l dx 3 



Here Hij is the transverse traceless gravitational wave 
term and A, Bi, Hl and Hrij are scalar degrees of 
freedom, two of which can be removed by gauge trans- 
formations. In Fourier space Bi = —kiB and i?r»j = 
(ktkj — Sij/3)HT- Often one uses longitudinal (or New- 
tonian) gauge with B = Ht = but we shall not use 
longitudinal gauge here. This is useful if we want to de- 
termine whether a given expression is gauge invariant. In 
a generic gauge, the g aug e invariant Bardeen potentials 
<f> and \I/ are given by [2l| 



f = A 
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-H 7 
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(A2) 
(A3) 



In longitudinal gauge they reduce to A = <3>, Hl = —<!>. 

From this we easily obtain the following expressions 
for the redshift perturbation [J [l(| , 
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1 + z 
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which leads to the density fluctuation in redshift space 
given in Eq. ([8]). To determine the volume perturbation 
we have to compute the derivative, 
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Here we made use of the fact that the Fourier transforms 
of Bi and Hxtj arc respectively 
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(A7) 



Using ^ = A + n l diX = X — d r X, we then obtain 
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The angular volume perturbation is, 
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To simplify the expressions we combine the terms 
(d 2 H T /dX 2 - 2dH T /dX)/k 2 of and 2/{k 2 r s )dH T /dX 
of 28r/r with Aoffr/(fcrs) 2 of the angular volume per- 
turbation to 
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We also use the gauge invariant velocity potential pj, |10|| 



V = v - -H T 
k 



( A8 ) so that 



(A10) 



Here An denotes the angular part of the Laplacian. The 
second integral in Eq. (|A8[) is the contribution from grav- 
itational waves which we shall not discuss further in this 
work. 

Putting it all together in Eq. (fl"8]l , using also Eq. ((29)) 
for the derivative of the perturbed redshift, we obtain the 
volume perturbation at fixed conformal time t or fixed 
background redshift z. However, we need to evaluate the 
volume fluctuation at a fixed, observed redshift which is 
related to the latter by 
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and the derivative of $ along the light ray, 
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With the help of these identities, the volume den- 
sity perturbation reduces to the gauge invariant expres- 
sion pop , where the gravitational wave contribution is 
omitted for simplicity. 



Appendix B: The contributions to the angular power spectrum 



In this appendix we express certain contributions to the total C' e s which are of particular interest for the discussion 
given in the text. We use the transfer functions for the concordance model given in eqs. (|45j) to (|48)l . 
Density fluctuation: 

Let us first consider the density term. The term of Tu in (|46p proportional to k 2 T^ largely dominates the integral. 
Its contribution is 

This integral only converges since T(k) decays like 1/k 2 for k > fc eq , where k cq is the (comoving) horizon scale at 
equal matter and radiation, see e.g. This integral is always dominated by the fluctuations on this scale, even at 
low I < £ cq (z) ~ 7rfc oq r(z s ). 
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Redshift- space distortion: 

The term T v (k) (j'/(krs)) 2 coming from d r (V ■ n) is the redshift space distortion. Since it is multiplied by k/H its 
dominant contribution behaves like the density term and is of the same order 
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Cross-term density-redshift space distortion: 
Also this term is of the same order as the previous two and even dominates at low 
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Lensing: 

The lensing contribution is in principle also of the same order. But since it probes the power spectrum truly at 
k ~ l/r(zs), it is largely subdominant at low I if compared to the previous contributions. 
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In the last equality we have used Limber's approximation [l7| 
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Cross-term density-lensing: 
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The other terms are 
Velocity: 
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Cross-term density-velocity: 
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In the regimes investigated in this work the gravitational potential terms are always subdominant and we do not 
write them down explicitly. 
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